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Abstract. Cayley's hyperdeterminant is a homogeneous polynomial of degree 
4 in the 8 entries of a 2 X 2 X 2 array. It is the simplest (nonconstant) polynomial 
which is invariant under changes of basis in three directions. We use elementary 
facts about representations of the 3-dimensional simple Lie algebra sl2(C) to 
reduce the problem of finding the invariant polynomials for a 2 X 2 X 2 array 
to a combinatorial problem on the enumeration of 2 X 2 X 2 arrays with non- 
negative integer entries. We then apply results from linear algebra to obtain 
a new proof that Cayley's hyperdeterminant generates all the invariants. 



1. Introduction 

In his famous 1845 paper on the theory of linear transformations, which became 
the foundation of classical invariant theory, Cayley [4] introduced the concept of 
the hyperdeterminant of a multidimensional array. He explicitly calculated the 
hyperdeterminant for the simplest case, an array of size 2x2x2, which can be 
represented in two dimensions by its two frontal slices: 



(1) X = 



Xwi Xlll 



Definition 1. Cayley's hyperdeterminant is the following homogeneous poly- 
nomial of degree 4 in the 8 entries Xijk of the 2x2x2 array of equation fl}: 

(~l _ 2 2,2 2,2 2,2 2 

- 2(x oo£ooi£iio2;m + zooo^oiozioi^m + ^ooo^oii^ioo^m 

+ a;ooi a^oio^ioi a^no + ^ooi^oii^ioo^no + ^oio^on^ioo^ioi) 
+ 4(x oo2;oiia;ioi2;iio + ^ooi^oio^ioo^m)- 

This polynomial has an interesting combinatorial-geometric interpretation. The 
first four terms have coefficient 1, and the subscripts correspond to the vertices of 
diagonals of the cube (configurations of dimension 1). The next six terms have 
coefficient —2, and the subscripts correspond to squares in the cube (configura- 
tions of dimension 2). The last two terms have coefficient 4, and the subscripts 
correspond to tetrahedra in the cube (configurations of dimension 3). These three 
configurations are illustrated by the dashed lines in Figure [1] 

Cayley's hyperdeterminant C is the simplest (nonconstant) polynomial in the 
entries of the 2x2x2 array X of equation ([1} which is invariant under unimodular 
changes of basis along the three directions. To make this idea more precise, we 
regard X as an element of the tensor cube C 2 (g> C 2 <g> C 2 of the 2-dimensional 
complex vector space C 2 . The group 6X2 (C) of 2 x 2 matrices of determinant 1 
acts on C 2 by matrix-vector multiplication, and this gives a component-wise action 

l 
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Figure 1. Geometric configurations in Cayley's hyperdeterminant 



of the direct product SL 2 (C) x SL 2 (C) x SL 2 (C) on C 2 <g> C 2 (g) C 2 . This action 
extends to the algebra of polynomials in the entries of X, and C is the simplest 
polynomial which is fixed by every element of the direct product. 

Ordinary determinants of square matrices can be characterized by a similar in- 
variance property. Matrices U G SL m (<C) act on rectangular m x n matrices A by 
left multiplication: A >-» UA. The First Fundamental Theorem of Classical Invari- 
ant Theory states that there exist nonconstant invariant polynomials in the entries 
of A if and only if m < n, and every invariant is a polynomial in the determinants 
of the m x m submatrices obtained by choosing m columns of A; see Procesi 
§11.1.2]. If we combine the left action of U £ SL m (C) with the right action of 
V e SL n (C), so that A H> UAV, then invariants exist for SL m (C) x SL n (C) if and 
only if m = n, and every invariant is a polynomial in det(A). 

We now summarize the results of this paper. In Section [2] we recall some ele- 
mentary results in the representation theory of the 3-dimensional simple Lie algebra 
512(C). We explain how the 9-dimensional semisimple Lie algebra 

sl 2 (C) 3 = s[ 2 (C)©sl 2 (C)©sl 2 (C), 

acts on the 8-dimensional vector space 

M 2 , 2 , 2 (C) =C 2 (g)C 2 ®C 2 , 

the tensor cube of the natural representation of sl 2 (C). We describe, using what 
are essentially the power and product rules from elementary calculus, the action 
of sl 2 (C) 3 on the algebra of polynomials on M 2 2 2 (C). The invariant polynomials 
are those which are annihilated by all Lie algebra elements (equivalently, fixed by 
all Lie group elements). For each degree d, the homogeneous polynomials form a 
finite-dimensional representation of s! 2 (C) 3 , and a well-known theorem implies that 
this representation is the direct sum of irreducible representations. We express the 
invariant polynomials as the elements in the kernel of a linear differential operator 
which represents the action of the Lie algebra on homogeneous polynomials, and 
from this we represent the invariant polynomials as the nullspace of a matrix. The 
domain of this linear map has a monomial basis in bijcction with the set of all 
2x2x2 arrays with non-negative integer entries summing to d and equal sums 
over the parallel 2x2 slices in the three directions. This reduces the computation 
of invariants to elementary combinatorics and linear algebra. 

In Section [3] we present explicit calculations for degrees 2 and 4. In degree 2, 
the matrix has size 6x4 and rank 4, so there are no invariants. In degree 4, the 
matrix has size 24 x 12 and rank 11; Cayley's hyperdeterminant C is a basis for the 
nullspace. Considering the powers of C, it follows that the dimension of the space 
of invariants is > 1 in each degree d which is a multiple of 4. 
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In Section |4] we compute the dimensions of certain weight spaces in the represen- 
tation of sb(C) 3 on the homogeneous polynomials of degree d. This is equivalent to 
the enumeration of 2 x 2 x 2 arrays with non-negative integer entries and constraints 
on the entry sums over the parallel 2x2 slices in the three directions. 

In Scction[5jwe apply a result on subspaces, reminiscent of the inclusion-exclusion 
principle, to a commutative diagram of injective linear maps between weight spaces 
in representations of s^C) 3 . This provides a different proof that the space of 
invariant polynomials has dimension > 1 in each degree d which is a multiple of 
4. We then use the representation theory of s^C) 3 to prove that the algebra of 
invariants is a polynomial algebra and is generated by Cayley's hyperdeterminant 
in degree 4. Hence there are no new invariants in higher degrees. 

In Section [6] we consider invariant polynomials in the entries of an array of size 
m x ni x • • • x rik under the action of SL ni (C) x SL, l2 (C) x • • • x SL Uk (C). The 
corresponding combinatorial objects are fc-dimensional arrays with non-negative 
integer entries and equal sums over the parallel slices in the k directions. 

In Section [7] we briefly summarize recent applications of Cayley's hyperdetermi- 
nant and provide some suggestions for further research. 



2. Representations of sl 2 (C) 



In this section we recall some elementary results in the representation theory of 
Lie algebras. Standard references are Jacobson [15] . Humphreys [14], de Graaf [5J, 
Erdmann and Wildon [9] . For an introduction to Lie theory, by which is meant the 
relation between Lie groups and Lie algebras, see Stillwell [24]. For the connection 
with classical invariant theory, see Procesi [22] . 

The 3-dimensional simple Lie algebra s^C) consists of the 2x2 matrices of 
trace over C with the Lie bracket operation [A, B] = AB — BA. This operation 
satisfies anticommutativity and the Jacobi identity: 



[A, A] = 0, [[A, B],C} + [[£?, C],A] + [[C, A],B] 
The standard basis of s [2(C) consists of these three matrices: 



0. 



(2) 



H 



1 
-1 



E 



F 





1 



In its natural representation, s^C) acts by matrix- vector multiplication on the 
two-dimensional vector space C 2 with this standard basis: 





' 1 " 




' " 


Xq = 





, X\ = 


1 



Lemma 2. The action of sl^C) on C 2 is given by the following equations: 

H ■ Xq = Xq, H ■ X\ = —X\. E ■ Xq = 0, E ■ X\ — Xq, F ■ Xq — X\, F ■ X\ = 0. 

In particular, xq and X\ are eigenvectors for H and H ■ Xi = {—l) l Xi (i = 0, 1). 

Lemma 3. If we regard xq and x% as indeterminates, then we can express the 
action of sfefjC) on C 2 by partial differential operators as follows: 



H 



Xq 



d 
8xq 



d 
dxi 



E 



Xq 



d 

dxi ' 



F = xi 



d 
8xq 
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Definition 4. We identify a 2 x 2 x 2 array X — (x^k) with an element of the 
tensor cube M2,2,2(C) = C 2 ©C 2 ©C 2 . We identify the entries with simple tensors, 
Xijk — Xi © Xj © Xk (i,j,k — 0,1). (Strictly speaking, since we regard Xijk as 
a coordinate function on M 2 ^ 2 ^ 2 (C), we should use dual basis vectors and write 
Xijk = x* © x* © x^, , but this distinction will not be important for us.) 

Definition 5. The Lie group SL 2 (C) x SL 2 (C) x 6X2 (C) acts on the vector space 
-^2,2,2 (C); the action is defined on simple tensors and extended linearly: 

{X, Y,Z) ■ (u®v®w) = (X • u) © (T -v)®{Z ■ w). 

As usual, we linearize the group action by considering the action of the Lie algebra 
S [ 2 (C) 3 = s( 2 (C) © sl 2 (C) © sl 2 (C) on M 2 ,2,2(C) defined by this equation: 

(A, B, C) ■ (u © v © w) = (A ■ u) ® v <g> w + u <g> (B ■ v) © w + u © v © (C • w). 

Lemma 6. The 8- dimensional vector space A/2, 2.2(C) is an irreducible representa- 
tion of the 9- dimensional semisimple Lie algebra s^C) 3 . 

Proof. A representation of a semisimple Lie algebra is irreducible if and only if 
it is isomorphic to the tensor product of irreducible representations of its simple 
summands. See Proposition 1.1 of Neher, Savage and Senesi [3T]. □ 

Definition 7. We write Hg,Eg,F^ (I — 1,2,3) for the standard basis of the £-th 
copy of sla(C) in s^C) 3 ; see equation j2|). 



Lemma 8. 


The basis of $1%(C 


) 3 acts on the basis o/Af 2 ,2,2(C) 


as follows: 


H t 


■ %Qjk — ^Ojfc; 


H 2 


XiQk = XiQk, 


H 3 - 


XijQ 


— XijQ 1 


Hx 




H 2 


XiXk Xixki 




XijX 


XijX 1 


E x 


• xojk = 0, 


E 2 


xmk = 0, 


E 3 - 


XijO 


= 0, 


E x 




E 2 


XiXk — XiQk ) 


E 3 - 


XijX 


= %ij0 1 


Fx 


%0jk — •&ljk 1 


F 2 


XiQk — XiXk 5 


F 3 . 


XijQ 


— •"ijl 1 


Fx 


■ xxjk = 0, 


F 2 


Xixk = 0, 


F 3 . 


XijX 


= 0. 


Proof. This is a straightforward calculation. For example, for I 


= 1 we have 


Hx 


• x 0j k = Hx ■ (x ( 


B> Xj <i 


>■) x k ) = {Hx ■ x ) 


g> Xj € 


5 x fe 


= Xojk , 


Hx 


• xijk = Hx ■ (xi ( 


g) Xj <i 


5 x k ) = {Hx ■ xx) 


8> Xj € 


5 x fc 


= -Xxjk, 


Ex 


• x aj k = Ex ■ {x C 


S> Xj (§ 


) x k ) = {Ex ■ x ) i 


5 Xj © 


Xfe = 


= 0, 


Ex 


• xijk = Ex ■ {xx C 


S> Xj (§ 


) x k ) = {Ex ■ xx) d 


5 


Xk = 


= Xojk, 


Fx 


• X jk = Fx ■ (x $ 


5 Xj (g 


Xk) = {Fx ■ xo) (S 


) © 


Xk = 


' Xxjki 


Fx 


■ xxjk = Fx ■ {xx $ 


5 Xj <g 


Xk) = {Fx ■ xx) <£ 


) Xj © 


Xk = 


-- 0. 



The other cases are similar. □ 
Definition 9. We consider the polynomial algebra on M 2j2 , 2 {C): 
P = C[xooo,xoio,xioo,xiio,x oi,xoii,xioi,x m ]. 
Lemma 10. A basis of P over C consists of the monomials, 

n^ijk rr e 000 rr e 001 rr e 010 rr e 011 rr e 100 rr e Wl rr e 110 

x ijk ~ x 000 x 001 x 010 x 011 x 100 x 101 x 110 x lll ' 

i,j,fe=0,l 

where the exponents eijk are arbitrary non-negative integers. 
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Definition 11. The degree of a monomial is the sum of its exponents: 

i,j,fe=0,l 

We write Pd for the homogeneous subspace of P spanned by the monomials of 
degree d. We identify P\ with M 2j 2,2(C), so that a basis of Pi consists of the 
monomials of degree 1, namely x 00 o, Zooi, xqiq, x Q u, x iao, xioi, x iio, Xm- 

Lemma 12. There is a bijection between the monomials of degree d and the 2x2x2 
arrays E = (eijk) of non-negative integers summing to d. 

Lemma 13. The polynomial algebra P is graded by the degree: 

S = ($P d , PdPe^Pd+e- 

d>0 

Lemma 14. We have Pd = S Pi, the d-th symmetric power of Pi. The action of 
an element D £ s^C) 3 extends to all basis monomials of P by the derivation rule 
D ' if 9) = {D ' f)9 + f{D ■ g). It follows by induction that 

D ' X ijk = e ijk x ij J k {D ' x ijk), 

and hence that 

i,j,k i' J f .k' 

= E ^OO™ ' ' ' {^i' j'k' x i'j'k' (P ' x i'j'k')j ■ ■ ■ x iri 
i' J' ,k' 

In particular, D ■ P d C P d for all D £ s[ 2 (C) 3 . 

Lemma 15. For every d > 1, the subspace Pd is a finite- dimensional representation 
of 512(C) 3 , and is therefore isomorphic to a direct sum of irreducible representations. 

Lemma 16. For every non-negative integer n, there is (up to isomorphism) a 
unique irreducible representation o/s^C) with dimension n+1, denoted V(n). This 
representation is generated by a vector v n with H ■ v n — nv n . With respect to the 
basis v n -2i (i = 0, 1, . . . , n), the action of shiC) on V(n) is given by 

H ■ v n -2i = (n-2z)v„_2i, 

E ■ v n -2i = (n-i+l)v n -2i+2 (i = 1,2, ...,n), E ■ v n = 0, 
F ■ v n -2i = (i+l)v n -2i-2 {i = 0, 1, . . . ,n-l), F ■ v_„ = 0. 

Definition 17. In the irreducible representation V(n) of 5(2 (C), the basis vector 
Vn-2i is a weight vector (that is, i/-eigenvector) of weight n—2i. 

Lemma 18. An irreducible representation of 3(2 (C) 3 is isomorphic to the tensor 
product V(a) (g> V(b) ® V(c) for some non-negative integers a, b, c. 

Lemma 19. The polynomials invariant under the group SL2{C) 3 coincide with the 
polynomials annihilated by the Lie algebra s(2(C) 3 . 

Lemma 20. A polynomial f £ Pd is invariant if and only if D ■ f = for all 
D £ 512(C) 3 . Equivalently, the invariant polynomials correspond to the summands 
of Pd isomorphic to V(0)®V(0)(g)V(0). Therefore, a polynomial f £ Pd is invariant 
if and only if Hi ■ f = and Eg ■ f = for 1= 1,2,3. 
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Lemma 21. The basis monomial 



t ijk ' 



is a simultaneous eigenvector for H 1 ,H 2 ,H 3 with eigenvalues 

Proof. For / = 1 we have 

Hi ■ x ijk = e-ijkx^' 1 (Hi ■ Xijk) = e iJ kx e l > k k ~ 1 (-iyx ijk = {-l) l e ijk x^ , 
and therefore 

//:• II •'•;,/ (E e °i k - E "•<) II 4* ■ 

The other two cases are similar. □ 

Definition 22. The weight space W(d; a,b, c) is the subspace of Pd spanned by 
the monomials which have eigenvalues a, b, c for Hi, H 2 , H 3 respectively. The 
zero weight space is W(d; 0,0,0). 



Lemma 23. The basis monomial 



II* 



e ijfc 
i-,3,k 



belongs to W(d; 0, 0, 0) if and only if 

ij'.fe j'.fe j.k i,k i,k i,j i,j 

That is, the 2x2x2 array (e^ ) of exponents satisfies the condition that in each 
of the three directions, the parallel 2x2 slices have equal sums. 

Lemma 24. If d is odd then the zero weight space W(d; 0, 0, 0) is the zero subspace. 
In particular, there are no invariant polynomials in odd degrees. 

Lemma 25. The actions of Ei, E2, E 3 induce these linear maps on weight spaces: 

Ei: W(d; 0,0,0) W(d;2,0,0), 

E 2 : W(d; 0, 0, 0) W(d; 0, 2, 0), 

£3 : W{d; 0, 0, 0) -> W(d; 0, 0, 2). 

Definition 26. We define a linear map 

£ d : W{d; 0, 0, 0) — ► W(d; 2, 0, 0) W(d; 0, 2, 0) © W(d; 0, 0, 2), 

by the equation £ d (f ) = (E 1 - f, E 2 ■ f, E 3 ■ f ) for all / e W(d; 0, 0, 0). 

Lemma 27. TTie invariant polynomials in Pd coincide with the kernel of £d- 

We can represent the linear map £d by the matrix {£ d] with respect to the ordered 
monomial bases of the weight spaces. The size [£d] is 

( dim W(d; 2, 0, 0) + dim W(d; 0, 2, 0) + dim W(d; 0, 0, 2) ) x dim W(d; 0, 0, 0). 

In fact the three dimensions in parentheses are equal; this follows by considering 
the automorphisms of s^C) 3 which permute the three summands. 
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3. CAYLEY'S HYPERDETERMINANT VIA LINEAR ALGEBRA 



In this section we show by direct calculation that every (nonconstant) invariant 
polynomial of degree < 4 is a scalar multiple of Cayley's hyperdeterminant. 

Wc identify monomials with sequences of exponents lexicographically ordered by 
their triples of subscripts: 



n 



„ e ijk 

Hjk 



|eooo, e ooi, eoio, eon, eioo, eioi, euo, em 



Within each weight space, we order the basis monomials lexicographically. 

Lemma 28. There are no invariant polynomials in degree 2. 

Proof. A basis of the zero weight space W(2; 0, 0, 0) consists of four monomials, 

00011000 00100100 01000010 10000001. 

which label the columns of the matrix [82]- Each nonzero weight space W{2\ 0, 0, 2) 
W(2; 0, 2, 0), W(2; 2, 0, 0) has a basis of two monomials which label the rows of [£2] 



01100000 




- 


1 


1 





10010000 




1 








1 


01001000 


m = 


1 





1 





10000100 





1 





1 


00101000 




1 


1 








10000010 




. 





1 


1 



The matrix has full rank, and so its nullspace is {0}. 



□ 



Theorem 29. In degree 4, the space of invariant polynomials has dimension 1; 
every invariant is a scalar multiple of Cayley's hyperdeterminant C . 

Proof. A basis of the zero weight space W(A; 0,0,0) consists of the 12 monomials 
in Figured Each nonzero weight space W(4; 0, 0, 2), W(A; 0, 2, 0), W(4; 2, 0, 0) has 
a basis of 8 monomials; see Figure[3J which also displays the 24 x 12 matrix [£4] (we 
use dot for zero). Figure 0] gives the row canonical form of [£4] (we omit zero rows). 
The rank is 11, and Cayley's hyperdeterminant is a basis of the nullspace. □ 

Corollary 30. The dimension of the space of invariant polynomials is at least 1 
in each degree d congruent to modulo 4- 

Proof. The existence of Cayley's hyperdeterminant C in degree 4 implies that there 
is at least one invariant polynomial C e in each degree d = 4e. □ 

4. Dimension formulas for weight spaces 

Our next goal is to prove that there are no new invariants in higher degrees; 
in other words, that every invariant is a polynomial in C . To do this, we need to 
prove that the lower bound of Corollary [30] is also an upper bound. The first step 
is to obtain dimension formulas for certain weight spaces in the representation of 
SI2 (C) 3 on the space Pd of homogeneous polynomials of degree d. 
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00022000 00111100 00200200 01011010 01100110 01101001 
02000020 10010110 10011001 10100101 11000011 20000002 

Figure 2. Monomial basis for zero weight space W(4; 0, 0, 0) 



01111000 
01200100 
02100010 
10021000 
10110100 
11010010 
11100001 
20010001 

01012000 
01101100 
02001010 
10011100 
10100200 
11000110 
11001001 
20000101 

00112000 
00201100 
01101010 
10011010 
10100110 
10101001 
11000020 
20000011 



.1.1.1 

2 . 1 
. 1 



2 1 
1 



Figure 3. The matrix [£ 4 ] 



1 -1 

.1 2 

. . 1 -1 

. . .1 2 

. . . .1 2 

1 -4 

1 .... -1 

1 ... -4 

1 . . 2 

1 . 2 



1 2 

Figure 4. The row canonical form of [£4] 
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Theorem 31. The dimension of the zero weight subspace W(d; 0, 0, 0) equals 

(000-0) -j- (d + 4) 2 {d 2 + 8d + 24) if d = (mod 4 ), 

384 

(000-2) — (d + 2) (d + 6) (d 2 + 8d + 28) i/ d = 2 (Vnod ^ 

384 

The dimensions of W(d; 2, 0, 0), W(d; 0,2,0) and W(d; 0, 0, 2) equal 

( 2 °0-0) ^7 d ( d + 4 ) 2 (d + 8) ifd = 0(mod4), 

384 

(200-2) ;Jr7 ( d + 2 ) (d + 6) (d 2 + 8d + 4) i/ d = 2 (mod 4 ). 

384 

T/ie dimensions of W(d; 2, 2, 0), VF(d;2,0,2) and VF(d; 0, 2, 2) eguaZ 

(220-0) — d(d + 4)(d 2 + 12d + 8) if d = (mod 4), 

384 

(220-2) (d + 2) (d 3 + 14d 2 + 28d - 24) if d = 2 (Wd ^ 

384 

The dimension ofW(d;2,2,2) equals 

(222-0) — d(d 3 + 16d 2 + 32d + 32) i/ d = (mod 4 ), 

384 

(222-2) — (d + 2) (d 3 + 14d 2 + 4d + 24) «/ d = 2 fmod ^ 

384 

In all cases, the dimension is if d is odd. 

Given non- negative integers d (the degree) and a,b,c (the weights), we consider 
2x2x2 arrays E = (e^fe) (i,j,k € {0,1}) of non-negative integer exponents 
satisfying the following equations: 

(D) eooo + eooi + eoio + eon + eioo + eioi + eno + em = d, 

(Ml) (e oo + eooi + e io + eon) - (eioo + eioi + eno + em) = a, 

(M2) (e oo + eooi + eioo + eioi) - (eoio + eon + e no + em) = b, 

(M3) (e oo + eoio + eioo + eno) - (eooi + eon + eioi + em) = c. 

These equations hold if and only if the corresponding monomial belongs to the 
weight space W(d; a, b, c); that is, the number of arrays E satisfying equations (|Dj) - 
(|M3p equals the dimension of W(d;a,b,c). Theorem |3"T1 gives formulas for these 
dimensions for certain values of a, b, c. These formulas are polynomials of degree 4, 
as expected since we have eight exponents and four constraints. 

Lemma 32. Consider 2x2 matrices (e^) with non-negative integer entries and 
specified row sums r , r*i and column sums Cq, Ci satisfying r + r% — c + Ci: 



eoo eoi 
eio en 



eoo + e i=7- , e o + eio = c , 

eio + en = ri, e i+eu=ci. 



(3) 

The number of such matrices equals min(r , fi, Cq, Ci) + 1. 



Proof. We have four variables and four constraints, but one dependence relation 
among the constraints, so we expect a 1-dimensional solution set. Without loss of 
generality, we can interchange the rows (resp. columns) and assume that ro < n 
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(resp. cq < ci); we can also transpose the matrix and assume that ro < cq. It is 



clear that since c\ — ro > 


ci - 


c > 


we have the particular solution 




eoo 


eoi 




' r 






eio 


en 




Co ci - r 





If u is any integer then we can preserve the constraints by adding u to the diagonal 
entries and subtracting u from the off-diagonal entries: 



eoo 


eoi 




u ro — u 


eio 


en 




c - u Ci — r + u 



This is another solution if and only if < u < ro. Hence the number of solutions 
is r + 1 = min(r , 7*1, Co, 0%) + 1. □ 

Lemma 33. For any integer k > 1 we have 



k k 

i=l j=l 



(i,j)) 2 = ^k(k+l)(k 2 + k + l). 



Proof. By induction on fc; the result is clear for k = 1. We have 



fc+l k+l 



XX (min(i, j))' 
i=i j=i 



k k 

XX 1 

i=l j=l 
1 



fc+l 

X 

i=l 
1 



( min(j, fc+l)) 2 + ( min(fc+l, j)Y 
j'=l 



-fc(fc+l)(fc 2 +fc+l) + -(fc+l)(fc+2)(2fc+3) + -fc(fc+l)(2fc+l) 



G 



6 



= -(fc + l)(fc + 2)(fc 2 +3fc + 3), 



using the formula for the sum of the squares from 1 to fc+l. 



□ 



We now come to the proof of Theorem [3TJ We prove the first two equations 
(1000-OP and (|000-2p : the proofs of the others are similar but slightly more compli- 
cated, and the details are not particularly enlightening. 



Proof. Equations (|M1|) (|M3|) imply that d is even, since ifa = fe = c = then each 
of the sums in parentheses equals d/2. Hence we assume that d = 2m. 

In equations (|Dl) - (|M3[) we write w,x,y,z for the row and column sums of the 
2x2 slice (eo^fe) with i — 0. We then have w + x = m and y + z = m, and 



eooo + eooi = w, 
eoio + eon = x, 
eooo + e io = 
eooi + eon = z, 



eioo + eioi =m — w, 
eno + em =m-x, 
eioo + eno =i rn — y, 
eioi + em = m - z. 



Suppose that w < x and y < z. Lemma [32] shows that 

• the number of 2 x 2 slices (eojk) is min(u;, y) + 1, and 

• the number of 2 x 2 slices (eyfc) is min(m— x, m—z) + 1 = min(w, y) + 1. 



CAYLEY'S HYPERDETERMINANT 



11 



Hence the number of 2 x 2 x 2 arrays is (min(w, y) + l) 2 . Any solution with w < x 
has a corresponding solution with w > x obtained by interchanging the slices (e^ofe) 
and (e^ifc)- Any solution with y < z has a corresponding solution with y > z 
obtained by interchanging (eyo) and (eyi). 

We first prove equation (|000-2|) : the case d = 2 (mod 4). We have m = 2k— 1 
where k = (d+2)/4. Since m is odd, we cannot have either w = x or y = z; hence 
all solutions are doubly paired. Thus the number of solutions is four times the 
number with w < x and y < z, and for this we apply Lemma [3"31 

fe-1 fc-1 fe-1 fe-1 

4 H Z! ( min ( w > rf + if^EI min(w+l, y+l) 2 

w—0 y—0 w—0 y—0 

= \k(k + 1)(A: 2 + fc + 1) = _L(d + 2)(d + 6)(d 2 + 8d + 28). 
3 384 

We next prove equation (|000-0j) : the case d = (mod 4). We have m — 2k where 

k = d/4. In this case we must also consider w = x and y = z, so we add 

fe-i fe-i 
2 ^ min(w+l,fc+l) 2 + 2^mm(k+l,y+l) 2 + min(fc+l, fc+1) 2 

w—0 y—0 

= \k{k + \){2k + 1) + (k + l) 2 = \{k + l)(4fc 2 + 5fc + 3), 

to the previous result, obtaining 

lk(k + l)(k 2 + k + l) + l(k + l)(4k 2 + 5fc + 3) = -i-(rf + 4) 2 (d 2 +8d + 24). 
3 3 384 

This completes the proof. □ 
5. Inclusion-exclusion for subspaces 

We recall a familiar formula from elementary linear algebra. If Ui and U2 are 
finite-dimensional subspaces of a vector space then 

(4) dim( Ui + U 2 ) = dim( U x ) + dim( U 2 ) - dim( U x n U 2 ). 

The next result generalizes equation (@} to an arbitrary finite number of subspaces, 
and is similar to the combinatorial formula for inclusion-exclusion on finite sets. 

Lemma 34. If U x , . . . , U n are finite- dimensional subspaces of a vector space then 

n n 

dim ( ^ ) < ( - !) r+1 £ dim ( ^1 n • • • n c/ ir ), 

i=l r=l l<ii<---<v<n 

where the inner sum on the right is over all (™) subsets {i x , . . . , i r } C {1, . . . , n}. 

Proof. The statement is false if the inequality is replaced by an equality: consider 
three distinct lines through the origin in the plane. The proof is by induction on 
n. The statement is clear for n < 2. We assume the statement for n and prove it 
for Using equation ((U) we obtain 



lim ( 



n+1 n 

Y, u *) = dim ((J2 Ui ) + f/ "+! 

»=1 i=l 



11 n 

dim ( £ 17<) + dim([/„+i) - dim ( ( ^ n C/„+i 

i=l »=1 
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Observing that 

n n 

dim((^f/i) nc/„+i) > dim( ^(C^nt^+i)), 

i=l i=l 

we obtain 

n+l n n 

dim(^f/ i ) < dim(^f/ i ) + dim(t/ n+1 )- dim ( £/; n t/ n+ i ) ) . 

i— 1 i—1 i—l 

We apply the inductive hypothesis to the last two sums, obtaining 

n 

5Z(-i) r+1 dim{U ii n---nU ir ) + dim(U n+1 ) 

r—1 \<ii<---<i r <n 
n 

-J2^ r+1 £ dim((f/ il nc/„ +1 )n---n([/ ir n?7„ +1 )). 

r=l l<ii <---<i r <n 

We separate the r — 1 terms of the first double sum, and simplify the second double 
sum using familiar properties of intersections: 

n 

]T dim(C/ j ) + I^(-l) r+1 Yl dim(U ii n---nU ir )+dim(U n+1 ) 

l<i<n r=2 l<i±<- ■ ■ <i r <n 

n 

_j2{-i) r+1 J2 dim(u il n---nu ir nu n+1 ). 

r=l l<ii<-..<i r <7l 

A slight rearrangement gives 

n 

dim((7 l )+E(- 1 ) r+1 E dim(U li n---nU lr ) 

l<i<n+l r=2 l<ii<-<v<n 

n 

+ ^ ( _ 1) (r+i)+i dim((c/ 41 n---nc/ J Jnc/„ +1 ). 

r=l l<ii < - <i r <n 

The first (resp. second) double sum corresponds to the subsets of size r (resp. size 
r+1) of the set {1, . . . , n+l} which exclude (resp. include) n+l, so we obtain 

n+l 

5>ir +1 ]T dim^n-.-ntU, 

r=l l<ii<---<i r <n+l 

and this completes the proof. □ 

We now consider a reformulation of this problem, in which we have a positive 
integer n and a collection of 2™ finite-dimensional vector spaces, 

{ Vi u i 2t ...,i n | < h,l2,- ■ -,i n < 1 }, 

corresponding to the vertices of an n-dimensional cube. We also have n2 n ~ 1 injec- 
tivc linear maps corresponding to the edges of the cube, 

where the hat indicates omission and the values of the indices are 
1 < k < n, (h, . . . ,i k , . . . ,i n ) e {O,!}"- 1 . 
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W(d; 0,2,0) 



F 2 



W{d; 0,0,0) ■ 



W(rf;0,2,2) 

F s 









F 2 



W(d;0,0,2) 

f 3 

f[ 




Ft 



Fo 



W(d;2,0,0) 



W(d, 2, 2, 2) 

F 2 



W(d;2,0,2) 



Figure 5. Linear maps among weight spaces in degree d 



Given any two of these vector spaces, we assume that all compositions of linear 
maps between the spaces give the same result; that is, the diagram is commutative. 

We can therefore identify each space Vi lt i 3t ,,, t i n with its image in Vb.o o, and so 

all of the spaces ^ii,i 2 ,...,j„ can be identified with subspaces of Vo,o,...,o- 

We define n vector spaces U\,...,U n by starting at the vertex (0, . . . , 0) of the 
rt-dimcnsional cube and following the n edges to the vertices 



Ui = v 



(1 < i < n), 



in which the subscripts on the right are except for 1 in position i. Given any 
r-element subset {zi, . . . , i r } C {1, . . . , n}, we write . . . , i r ) for the element of 
{0, l} n which has 1 in positions i\, . . . , i r and elsewhere. Our assumptions allow 
us to make the following identifications: 



u i± n • ■ • n Ui, 



V x(h,...,ir)- 
Lemma [34] then implies that 

dim(im(/W, )+...+im(/(;l )) < 

n 

£(-1)^ E dim(y x(il ,..., ir) ). 



(5) 



l<?i < ■ ••<i r < 1 n 



Example 35. We consider n = 3 and identify the 8 vertices of the cube with the 
following weight spaces in degree d defined in Section [4] 

W(d; 0,0,0), W(d;2,0,0), W(d; 0,2,0), W(d;0,0,2), 
W(d;2,2,0), W(d;2,0,2), W(d; 0,2,2), W(d;2,2,2). 

The representation theory of sl2(C) shows that the action of the basis elements 
^1,^2,^3 on the homogeneous polynomials of degree d gives injective linear maps 
between these weight spaces as illustrated in Figure [SJ The invariant polynomi- 
als are the nonzero elements in the irreducible summands ^(0) <g) V(0) g) V(0), 
and the number of these summands equals the codimension, in the zero weight 
space W(d; 0,0,0), of the sum of the images of the weight spaces W(d;2, 0,0), 
W(d; 0, 2, 0), W(d; 0, 0, 2) under the actions of Fx, F 2 , F 3 respectively. That is, 
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• We start with the entire zero weight zero space W(d; 0, 0, 0). 

• We factor out the images of vectors of weight (2,0,0) or (0,2,0) or (0,0,2) 
by the action of F-y or F2 or F3 . 

• The vectors that come from weight (2,2,0) or (2,0,2) or (0,2,2) by the action 
of F\, F2 or Fi, F3 or F2, F3 have then been factored out twice, so we must 
add those dimensions back in. 

• But then the vectors that come from weight (2,2,2) by the action of F\ , F2 , F3 
must be factored out again. 

The dimension formulas from Section [4] with equation ([5]) give 

dim W(d; 0, 0, 0) - dim W(d; 2, 0, 0) - dim W(d; 0, 2, 0) - dim W(d; 0, 0, 2) 
+ dim W(d; 2, 2, 0) + dim W(d; 2, 0, 2) + dim W(d; 0, 2, 2) - dim W(d; 2, 2, 2) 




1 if n = (mod 4) 
otherwise. 



Combining this with Lemma [MJ this gives another proof of Corollary [3DJ the 
dimension of the space of invariants is > 1 in degrees d = (mod 4). 

Theorem 36. Every polynomial in the entries Xij k of the 2x2x2 array X = (xij k ) 
k = 0, 1), which is invariant under changes of basis with determinant 1 along 
all the three directions, is a polynomial in Cayley 's hyper determinant. 

Proof. It remains to use the representation theory of Lie algebras to show that 
inequality ([5]) becomes in fact an equality in the situation of Example [35J We know 
that the space Pd of homogeneous polynomials of degree d is completely reducible 
as a representation of the semisimple Lie algebra 5(2 (C) 3 , and that the irreducible 
summands are tensor products V(a) ® V(b) (E> V{c) of irreducible representations 
of 512(C). Since the weight spaces in the tensor factors have dimension 1 as rep- 
resentations of sl2(C), it follows that the weight spaces in the tensor product have 
dimension 1 as representations of 3(2 (C) 3 . Inequality ([5]) is obviously an equality 
when all the dimensions are 1, and this completes the proof. □ 

6. General multidimensional arrays 

We consider a fc-dimensional array of size n\ x 71,2 x • • • x n k : 

X = (x ili2 ... ik ) (l<k<ni, 1 < i 2 < n 2 , 1 < ik < n k ). 

(The smallest index is now 1, not 0.) We consider an extension of determinants to 
these arrays, using a combinatorial approach based on the representation theory of 
the special linear Lie algebra sl„(C). As usual we write C" 1 , C™ 2 , . . . , C nk for the 
complex vector spaces with dimensions m, n 2 , ■ • ■ , nk and standard bases 

e-^ (h = 1, ...,ni), (i 2 = 1, . . . ,n 2 ), efj (i k = 1, . . . , n k ). 

A tensor of order k is an element of the tensor product 

C ni,tia,...,t»» _ C m ^ C n 2 (g, . . . g, . 

Lemma 37. Every element of C™ 1 '™ 2 ' '"'™ fc is a finite sum of elements of the form 

vi®v 2 ®---®v k («! eC ni ,v 2 eC n2 ,...,v k eC"'). 

A basis for C™ 1 '™ 2 ' '"'" fc over C consists of the n\U2 ■ ■ ■ n k simple tensors 
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Every tensor of order k can be expressed uniquely in the form 

ni 112 rife 
ii «2 ifc 

A fc-dimensional array consists of the coefficients of a tensor of order k with 
respect to the basis of simple tensors: 

x = ( x ii,i2,...,ik) (ii = 1, ...,ni; i 2 = 1, ...,n 2 ; ...; ifc = l,...,rife). 

If Mi , M 2 , • • • , Mfc are linear operators on C" 1 , C™ 2 , . . . C nk then, with respect to 
the standard bases, we identify Mi with an ni x ni matrix for I = 1, 2, . . . , k: 

M e = (m ( i f) (mg ) eC;i,j = l,...,n / ). 

The action of a fc-tuplc of operators M = (Mi, M2, . . . , Mfc) on a simple tensor in 
C" 1 '" 2 '— is given by the equation 

(6) (Mi, M 2 , . . . , M fe ) ■ (v! ® w 2 ® ■ • ■ ® v k ) = Miwi ® M 2 v 2 <8> • • • ® M k v k . 

We introduce nin 2 • • -n k indeterminates corresponding to the entries of X: 

£ii,i 2 , ...,»* (ii = 1, . . . , ni; 12 = 1, • • • , n 2 ; . . . ; i fe = 1, . . . , n fc ). 

We consider the polynomial algebra in these indeterminates over C: 

C[xi u i 2 ,...,i k I «i = 1, • • ■ ,ni; i 2 = 1, • • ■ ,n 2 ; . . . ; i k = 1, . . . ,n fe ]. 

For £ = 1, 2, . . . , fc the action of M( on an indeterminate corresponds to its action 
on the standard basis vectors in C™' : 



=J2 m if e i e) M C X h,-,h,-,jk=J2 m ijl Xi ^ 
i=l i=l 



From this we obtain the action of M = (Mi, M 2 , . . . , M k ) on an indeterminate: 

ni n 2 n k 

(Mi, m 2 , . . . , M fc ) . x hj2 ,... ijk = E E • • • E m Si m S2 • • • 

ii=li2=l »fe=l 

This action of A'/ = (Mi, M 2 , . . . , Mfe) extends to an action on polynomials: 

M ■ /(aJll...l, • • • , %jij 2 ...j k j • • • 7 x n 1 n 2 ...n k ) = 

/(M • X11...1, . . . , M • Xj 1 j 2 ...j k , . . . , M 
Definition 38. The polynomial / £ C[xj li i 2i ... i i fc ] is invariant if 

det(M/) = 1 (^ = 1, . . . , fc) => M • / = /, M = (Mi, M 2 , . . . , M fe ). 

The n x n complex matrices of determinant 1, with the usual operation of matrix 
multiplication, form the special linear group SL n (C). Finite-dimensional represen- 
tations of SL n (C) can be studied in terms of the Lie algebra sI„(C), which consists 
of all n x n complex matrices of trace 0; the bilinear product is the Lie bracket 
[A, B] = AB - BA. The standard basis of s[„(C) consists of 

• the matrix units Uij for i ^ j with entry 1 and other entries 0, 

• the diagonal matrices Hi — Uij — CA+1^+1 for z = 1,2,..., n— 1. 

The simple root vectors are the matrix units Ei = t/^i+i for i = 1, 2, . . . , n— 1. The 
natural representation of st„(C) is its action on C" by matrix- vector multiplication. 
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Lemma 39. In the natural representation of sln(C) we have 



ej-i if j = i+l 
otherwise. 



if j = i 

Hi ■ ej = ^ —ej if j = i+1 Ei ■ e 3 = 

otherwise, 

We consider the action of the semisimple Lie algebra 

k 

(8) 0si n ,(c)= 5 [ ril (C)e5i„ 2 (c)©---esU(c) ) 

on its irreducible representation C™ 1 ' ™ 2 '' ■■• Tlk , the tensor product of the natural 
representations of its simple summands. For £ = 1, 2, . . . , k we write H\ , E\ ' for 
the elements Hi, Ei 6 s\. ne {C). Combining equations ((6]) and (0) with Lemma [39] 
we obtain the action of H\ and £^ on the indeterminates Xj 1 j 2 ...j k . 

Lemma 40. For £ = 1,2, ... ,k and i = 1,2, . . . , nt—1 we have 

*3i,h,...,h ifjt=i 
Hi ' x ji,j2,...,j k = ) ~ x jt,h,—3h if ' 3t = 

otherwise, 

F (<) _ _ I x 3i,h,-,h-i,-j k ifje = i+l 

a i ' Jj jl-j2----,jk ) „ 

I (J otherwise. 

The action of a Lie algebra L on a tensor product V ® W of representations is 
given by the derivation rule: 

x ■ (v <S> w) = [x ■ v) <8> w + v <g> (x ■ w) (x S L, v S V, w € W). 

We identify the <i-th symmetric power iS d V of the representation V with the space 
of homogeneous polynomials of degree d on a basis of V. It follows by induction on 
d that the action of L on S d V is given by the following equation: 

p 

x ■ «uf 2 ■ ■ ■ V f>) = ^vl 1 ■ ■ ■ (x ■ uf) • • ■ 

i=l 

P P 

i=l i=l 
We apply this to 

k ni ri2 n k 

£=i ii=ij2=i jfc=i 

Some equations will be clearer if we write a monomial as follows: 



m ii2 «*: 

TT TT TT T, e 3u'2 -ifc _ „ei...i ._, e ii--ife . „ e "i - "fc 

1111 11^152...^ I ... *. X »l---«fe 

ii 32 jk 
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Lemma 41. For £ = 1, 2, ... ,k and i = 1, 2, . . . , n^ — 1 we have 

ttW . / ei...i e n - jfc . . . e "i--"fc\ _ 

■"» V x l-..l ji-.-jk x rn...n k ) — 

"' / j \"je-i ~ "je-i+l) e ji---j<!---jk x l.'.'.l " ' X ji...j e ...j k ' " x n 1 ...n k , 

ii=i ifc=i 

u h---3k 

«l n k 

E...ST X. ., , p . . . x ei...i . . . ~ e H-U- l -ih +1 , . . -Vi-it-ik' 1 . . . _ e »i-»fc 
2^ e Jl---Jf---Jfc X l...l ^ji-je-l-jk ^ji-je-jk x n 1 ...n k , 

ii=i ifc=i 



V x l...l • L h—3k x n 1 ...n k ) — 



where 5ij is the Kronecker delta (5a = 1, <5.y = for i =/= j). 

Lemma 42. For every I = 1, 2, . . . , k and i = 1, 2, . . . , ng—1, the monomial 



J l...l ■*'TH...f»fc 5 



is an eigenvector for H\ with eigenvalue 

ni _ — _ "i _ — _ 

" X] e h---i---3k x " ■ x " ■ x e Ji---*+i---jfe' 

ji=i ifc=i j'i=i ifc=i 

w/iere i/ie hat denotes omission. 

The space of homogeneous polynomials of degree d has the basis 



"l n k 



Ji = l Jfe = l 

Definition 43. A monomial x^ 1 ^ 1 ■ ■ ■ x^.'.'.n* has weight zero if it has eigenvalue 
for every H\ with £ = 1,2, ... ,k and i = 1,2,..., ng—1; that is, 

"l - — - »fc "i - — - . "fc 

^ ' ' ' X/ ' ' ' X/ e h-i-3k = E ' ' ' E ' ' ' X/ e ji-»+l-jfc- 

ii=i jfe=i ii=i je j k =i 

The zero weight space of degree d consists of the monomials of weight zero. 

Definition 44. Let E = (e^...^) be an array of size n\ x ri2 x • • • x with non- 
negative integer entries. A slice of E is a (k— l)-dimcnsional subarray obtained by 
fixing one subscript; for every £ = 1,2, ... ,k we can set ie = 1, 2, . . . , ng and obtain 
ng slices of size n\ x • • • ng ■ ■ ■ x n^. We call E an equal parallel slice (EPS) 
array if for every £ = 1, 2, . . . , k the ng slices in direction £ have the same entry 
sum. That is, for each i the following sum does not depend on j: 



ii=l if ifc = l 

Lemma 45. ^4 &oszs for the zero weight space in degree d consists of the monomials 
whose arrays of exponents are EPS arrays. 

We write W(d; a\, . . . , a„_i) for the vector space with basis consisting of the 
monomials with degree d and eigenvalues (ai, . . . , a„_i) for Hi, ... , H n -i as in 
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Lemma |42] In s(„(C) the brackets of Hi and Ej are given by the formulas 

!2Ej i£i = j 
-Ej if j = i-loij = i + l 
otherwise. 

It follows that the actions of E\, . . . , E n _i induce the following linear maps: 

Ex: W(d;0,... ,0) —>W(d; 2, -1,0,. ..,0,0), 
W{d; 0, . . . , 0) — > W(d; -1,2,-1,..., 0,0), 
W(d;0,...,0)^W(d;0,-l,2,...,0,0), 



E 2 
E, 



E n -i : W(d- 0, . . . , 0) — > W(d; 0, 0, 0, ... , -1, 2). 
The weights appearing on the right are the rows of the Killing-Cartan matrix, 

f 2 if * = j 

iff™- 1 ) = (/ty), «y = < -1 if j = * - 1 or j = i + 1 

[ otherwise. 

We write w^ 1 , . . . , U!^!^ 1 for the rows of K^ n ~^ and form the linear map 

n-i 

£ = (Bi,...,£ n _i): W(d;Q,...,0) — )• ® W(d;«;f 

»=i 

We apply this to the semisimple Lie algebra ©. We first combine the spaces 
W(d; 0, . . . , 0) for each summand into the zero weight space of Definition (|43|) : 

"l-l n k — 1 

z = w^(d ; fv^7o) n • • • n W^(d;6T^~0). 

We then combine the linear maps E for each summand into the single linear map 

k ni — 1 

(9) £ = {E^\. . . , : Z — > ® ® W(d; w^). 

i=\ i=i 

Theorem 46. The invariant polynomials in degree d for the n\ X ■ ■ ■ X r%k array 
X = (a;i 1 ..-i fe ) are the (nonzero) elements of the kernel of the linear map ([9]). 

In degree d there are no monomials of weight zero unless d is a multiple of 
N = LCM(ni, . . . , 7ifc); hence invariants can only exist in degrees d = (mod N). 

We expect that the dimension of the zero weight space in degree d (equivalently, 
the number of EPS arrays with entry sum d) is a polynomial in d. Since we have 
ni ■ • • rife exponents, with one constraint on the degree and (ni— 1) + • • ■ + {n^ — 1) 
constraints on the parallel slices, we make the following conjecture. 

Conjecture 47. Let k and rii, ri2, ■ ■ ■ , Uk be positive integers. The dimension of 
the zero weight space in degree d is given by a family of polynomials of degree 

k k 



j| tij - y] ni + k - 1. 



i=i t=i 



CAYLEY'S HYPERDETERMINANT 



19 



7. Conclusion 

Modern interest in Cayley's hyperdeterminant and its generalizations was re- 
vived by the famous paper of Gelfand, Kapranov and Zelevinsky |10j ; see especially 
Proposition 1.9 on page 234. The same authors developed this subject in great 
depth, using the techniques of algebraic geometry, in their monograph 

A closely related topic, of great importance in applied numerical linear algebra, 
is the problem of computing the rank of a fc-dimensional array. When k = 2, this 
problem has an efficient solution using Gaussian elimination, but for k > 3 it has 
been shown by Hastad [12j to be NP-complete. A comprehensize survey on tensor 
rank and algorithms for tensor decomposition has been given recently by Kolda and 
Bader 16J. Cayley's hyperdeterminant was rediscovered in the 1970's by Kruskal 
[17], and is sometimes called Kruskal's polynomial by applied mathematicians; see 
ten Berge [55] and Martin [3U] for an explanation of how it can be used to compute 
the rank of a 2 x 2 x 2 array. Two recent related papers are de Silva and Lim [5] 
and Stcgcman and Comon [23 . 

Invariant polynomials on arrays of size 2 x 2 x • • • x 2 (k factors) have been studied 
by theoretical physicists working on quantum computing; see Luque and Thibon 
[18] [19] , Djokovic and Osterloh [8] . For the combinatorial-geometric aspects of this 
problem, see Huggins et al. [13] . These invariants can be regarded as noncommuta- 
tive analogues of classical invariant theory (for a survey see Dixmier [7]): the 19th 
century invariant theorists studied the irreducible representations V(d) = S V(l) 
of sl2(C), and replacing the symmetric power by the full tensor power gives the 
vector space of arrays of size 2 k . It is an open problem to extend the methods 
of the present paper to these arrays. It would be very useful to have a complete 
description of the structure of the space of homogeneous polynomials as a sum of 
irreducible representations of the semisimple Lie algebra; one possible approach to 
this problem has been developed by Adsul and Subrahmanyam pQ. 

The objects that we call equal parallel slice (EPS) arrays are examples of contin- 
gency tables, which are important in combinatorics and statistics. For asymptotic 
formulas for the enumeration of these objects, see Barvinok [3J. In closing, we 
mention the intriguing applications of Grobner bases and hyperdeterminants to 
mathematical genetics; see Allman and Rhodes [2], especially page 146. 
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